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Cl, Abstract 

We outline the application of the truncated conformal space approach (TCSA) to perturbations of 
SU(2)k Wess-Zumino-Witten theories. As examples of this methodology, we consider two distinct 

> 

■^j- perturbations of SU(2)\ and one of SU(2)2- SU(2)\ is first perturbed by its spin-1/2 field, a 

00 

model which is equivalent to the sine-Gordon model at a particular value of its coupling f3. The 

O 

sine-Gordon spectrum is correctly reproduced as well as the corresponding finite size corrections. 

o 

We next study SU (2)i with a marginal current-current perturbation. The TCSA results can be 

matched to perturbation theory within an appropriate treatment of the UV divergences. Finally, 

^ we consider SU(2)2 perturbed by its spin-1 field, which is equivalent to three decoupled massive 

H 

Majorana fermions. In this case as well the TCSA reproduces accurately the known spectrum. 
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I. INTRODUCTION 



Correlations in one dimensional systems typically require non-perturbative, non-mean 
field techniques to access the underlying physics. These techniques often trade on low energy, 
field theoretical reductions of the system. Examples include bosonization [1], conformal field 
theory [2], integrable field theory [3j, the Bethe ansatz jl], and the truncated conformal 
spectrum approach (TCSA) [5J. The last approach, unlike the aforementioned techniques, is 
able to deal in principle with any one dimensional field theory in an exact numerical manner. 
In that sense the TCSA is similar to the density matrix renormalization group (DMRG) [B], 
but the framework where it is formulated is field theoretical and not discrete quantum lattice 
systems. 

TCSA deals with models whose Hamiltonians can be represented in the form, 

H = H C FT + (1) 

where He ft is the Hamiltonian of some conformal field theory and $ is an arbitrary per- 
turbation. The approach employs the Hilbert space of the conformal field theory as a 
computational basis and exploits the ability to compute matrix elements of the perturbing 
field in this same basis using the constraints afforded by the conformal symmetry. It was first 
employed by Yurov and Zamolodchikov in studies of massive perturbations of the critical 
Ising model [7j and the scaling Yang-Lee model [5]. 

Since its introduction it has been used in a large number of instances and to some degree 
has become a standard tool. It has been used to study perturbations of the tri-critical 
Ising [8j [9], the 3-states Potts model [10J, bosonic (c = 1) compactified theories [11], the 
sine-Gordon model [121 02], and perturbations of boundary conformal field theories [T%] . 
Spectral flows between different conformal field theories were addressed in Ref. [TH], while 
the correctness of the thermodynamic Bethe ansatz equations was checked in Refs. [T6| IT7]. 
Finite-size corrections to the mass spectra have also been analyzed in Refs. [T8| [19] . By 
replacing the conformal field theory with an integrable field theory, the approach can also 
study models of perturbed integrable field theories. Matrix elements of the perturbing field 
are then computed in the form factor bootstrap approach [2D] . 

While the TCSA approach is extremely flexible in the models it can attack, it is in 
practice limited to perturbations of conformal field theories with small central charge (c ^ 1). 
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For theories with large central charge, the underlying conformal Hilbert space is large and 
becomes numerically burdensome to manipulate. The difficulty has recently been partially 
ameliorated with the development of a numerical renormalization group for the TCSA [2TL 
[25] . This renormalization group permits large Hilbert spaces to be dealt with piecewise 
making the numerics manageable. Using this renormalization group, the excitonic spectrum 
of semi-conducting carbon nanotubes was studied [23] (here the underlying conformal field 
theory had c = 4) as were large arrays of coupled quantum Ising chains [22] (here the 
underlying conformal field theory had c ~ 30 — 50). 

The TCSA approach, as designed, focuses on accurately computing the properties of 
the low energy states. However when combined with an NRG together with a sweeping 
algorithm not dissimilar to the finite volume algorithm of the DMRG [6], the TCSA can 
compute the properties of states over a wide range of energies. This was demonstrated in 
[23] where the level spacing statistics were studied in crossing over from an integrable to a 
non-integrable model. 

To the best of our knowledge the TCSA has not been applied previously to perturbations 
of Wess-Zumino-Witten (WZW) models. WZW models are non-linear sigma models whose 
field g lives on a group manifold G. They possess topological terms, the Wess-Zumino term, 
whose action is quantized with the consequence that its coupling constant k is constrained to 
be a positive integer. The consequence of the topological term is to make the sigma models 
conformal with the affine Lie algebra associated with G, spectrum generating. Affine Lie 
algebras typically have richer structure than the Virasoro algebra and are consequently more 
difficult to treat with the TCSA. Specialized code needs to be developed in order to treat 
such models. Here we report the development of such code for the study of perturbations 
of SU(2) k WZW models. 

Perturbed SU{2) k WZW models are interesting physically primarily because they are able 
to represent the low energy structure of spin chains [25]. The most important example here is 
the spin- 1/2 Heisenberg chain, whose low energy behavior is governed by 577(2) i perturbed 
by a marginally irrelevant current-current interaction. Adding different perturbations to 
SU (2)i leads to different variants of the Heisenberg model. The Hamiltonian 



H = H SU {2) 1 + 9 / dx J R - J L + h / <ix (01/2,1/2 01/2 -1/2 - 01/2-1/2 01/2,1/2) , (2) 



where 0i/2,±i/2 are the spin- 1/2 fields in S77(2)i, is the low energy reduction of the dimerized 
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J\ — J2 Heisenberg model 

H = Ji (! + <K-!) n ) ^ • 5 n+1 + J 2 S n S n+1 . (3) 

n 

The couplings of the two models are related via g oc J2 — J2C and h oc 5. This model 
has been studied intensely both field theoretically with RG analyses [27] and with DMRG 
[281 [29] . These methodologies do not currently agree on how the spin gap depends upon the 
dimerization parameter 5. It is one of the aims of our work to set up the framework under 
which disputed questions surrounding this model can be addressed. 

SU(2)k for k > 1 WZW theories are also of considerable interest as they are the low energy 
reductions of families of spin-fc/2 spin chains with finely tuned, local interactions [26]. They 
have also been shown more recently to represent Haldane-Shastry type spin chains [30] with 
longer range interactions [3"TH3"3"] . In both cases, it is of interest to understand how SU(2)k 
WZW behaves in the presence of relevant and marginal perturbations. And more generally, 
because SU(2)k WZW theories are multicritical with many possible relevant perturbations. 
Actual spin chains are likely to be realized only in the vicinity of these critical points rather 
than precisely at them. 

The outline of this paper is as follows. In Section II we present briefly the framework 
for the TCSA where perturbations of SU(2) k can be studied. In the next three sections 
we present applications of the TCSA to perturbed 577 (2)^ models, establishing how the 
methodology works. In Section III we examine the SU(2)i+Tr(g) model, which is equivalent 
to the sine-Gordon model. This allows us to compare our numerics with known analytic 
results. In Section IV we study the SU(2)% + Jl - Jr model, which corresponds to a marginal 
current-current perturbation of SU(2)\. We show here that one can isolate the UV divergent 
behavior, essential for extracting the universal behavior of spin chains described by SU(2)k 
with a marginal perturbation. Finally, in Section V, we consider SU(2) 2 perturbed by 
Tr(g) 2 . This provides a useful benchmark of our methodology as the theory is equivalent to 
three non-interacting massive Majorana fermions. We close in Section VI with conclusions 
and a discussion of future directions. 
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II. TCSA FOR SU(2) k WZW MODELS 



In this section we describe the application of the TCSA [S] to deformations of the SU (2)^ 
Wess-Zumino-Witten model [Ml 155]. We thus start by considering Hamiltonians of the form 



H = H su{2h +g [ dx$(x) 
Jo 



(4) 



Here $ is a spin singlet combination of (highest weight or current) fields of the WZW 
model. The Hamiltonian is defined on a circle of length R. Adding the time coordinate the 
underlying space-time is an infinite cylinder with circumference R. 

The first step in the TCSA is to characterize the unperturbed theory, Hsu(2) k - This 
theory provides the computational basis of the TCSA numerics. Hsu(2) k has central charge 
c = 3k/(k 4- 2) (k = 1, 2, . . .), and can be written a la Sugawara in terms of the SU(2)k 
currents: 

2tt 



H 



SU(2) k 



R 

2tt 
~R 



c 
12 



m 

-l(2 7° 7° 4- 7 + J~ 4- 7 _ 7 

1 \ m —m 1 "m"—m ' " m" —ml 



C 

12 



(5) 



The left-moving currents obey the algebra, 

km 



\J° J°l 

\J° J±] 
l u mi °n J 

\j~rni Jn ] 
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±J. 



n+m,0j 

± 



(6) 



2J m +n + kmd m+nfi , 



with the right moving currents J^f 1 obeying the same algebra. 

The field content of Hsu(2) k consists of k + 1 primary fields, <p 8 ,m=-8,-,si forming spin 
s = 0, • • • , k/2 representations. The conformal weight of the spin s primary field is given 
by A s = s(s + l)/(fc + 2). The primary fields and the WZW currents are the basic tools 
to construct the Hilbert space of the unperturbed theory, providing the first ingredient of 
the TCSA. The Hilbert space can be written as a tensor product of its holomorphic (left) 
and anti-holomorphic (right) degrees of freedom. For SU(2)}., the left and right sectors of 
the Hilbert space each have k + 1 modules, where each module is associated to one primary 
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field. The module consists of a highest weight state, \s,s) = </> S)S (0)|0), together with an 
infinite tower of descendant states: 



These states ([7]) are eigenstates of the Virasoro operator L , and the third component 
component of the current Jj: 



The quantity J2i n « i s called the Kac- Moody level, or simply the level, of the descendant 
state. 

The set of states (f7| is not yet a basis of the Hilbert space since it is over complete 
and contains null states. In order to form a complete orthonormal basis, we tackle each 
Kac-Moody module separately. We do so in an iterative fashion. At each step we have a 
set of non-zero norm states which are linearly independent (at the beginning this set will 
consist solely of the highest weight state). We next add a new descendant state to this set, 
compute the matrix of scalar products of states in the expanded set (the Gramm matrix), 
and find its determinant. If it is non-zero, the new state is added to the list; otherwise it is 
discarded. We then move to the next descendant in the tower of states in increasing order 
in the level. The process ends when all the states up to a given level have been considered. 
This procedure yields a complete set of states that we can easily orthonormalize to obtain 
a basis. In order to optimize this procedure we take into account the following properties: 

• states with different L quantum number are independent; 

• states having different spin (eigenvalue of J°) are independent; 

• we discard the states with null norm; 

• we act only with level currents, Jq' , directly on the highest weight state \s,s) in 
the module. 




7li = 0, 1, 2, • • • <2j = 0, ±. 



(7) 



L (T™ M . . . J%) \s, s)= A. + £ n, (JX U . . . JX) K *>; 



(8) 



6 





5Z7(2)i 




SU(2) 2 






level 


spin 


spin 1/2 


spin 


spin 1/2 


spin 1 





1 


2 


1 


2 


3 


1 


4 


4 


4 


8 


7 


2 


8 


10 


13 


20 


19 


3 


15 


18 


28 


46 


40 


4 


28 


32 


58 


94 


83 


5 


47 


52 


112 


178 


152 


6 


76 


86 


206 


324 


275 


7 


119 


132 


359 


564 


468 


8 


181 


202 


611 


948 


786 


9 


271 
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1272 


10 


397 
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1611 
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2026 


11 


571 


622 


2529 


3886 


3145 



TABLE I: The dimensions of the Verma modules of SU(2)\ and SU{2)2 at different levels. 

To demonstrate that the above method amounts to a numerically intensive task, we present 
in Table I the number of states per level for the two modules of S77(2)i and the three 
modules of SU (2)2. 

Once the chiral sector of SU(2)k has been obtained, the total Hilbert space is constructed 
as a tensor product of the isomorphic holomorphic and antiholomorphic sectors. These tensor 
products are diagonal in the modules, i.e. left moving spin s states are only tensored with 
their right moving spin s counterparts. In forming these tensor products, we group the 
states by their value of Lorentz spin, (L — L ), and ^-component of SU(2) spin, Jq + J°. 
Recall that the Lorentz spin is proportional to the momentum carried by the corresponding 
state. In Table [H| we present the number of cumulative states up to a given level in S77(2)i 
and S77(2) 2 with vanishing Lorentz spin and 577(2) spin. 

Once the computational basis has been constructed, the TCSA requires the evaluation 
of matrix elements of the perturbing operator in that basis. For this purpose one uses the 
commutation relations of the fields, s , m (O), with the current modes, Jq j± , which are given 
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SU(2) 1 level 
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SU(2) spin 


SU(2) 2 level 
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75 


25 
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24 
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120 
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1309 


381 
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6839 


1595 


11 


133123 


32021 
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185111 


38665 



TABLE II: The cumulative dimensions for non chiral spaces at given level with total Lorentz spin 
(mom.) zero (first column) and both Lorentz spin and SU(2) spin S z zero (second column). 



SU(2) 1 




Co,0;l/2,-l/2;l/2,l/2 1 


Cl/2,l/2;l/2,l/2;0,0 1 


SU{2) 2 




Co,0;l/2,-l/2;l/2,l/2 1 
Co,0;l,-l;l,l 1 
Cl/2,l/2;l/2,l/2;0,0 1 
Cl/2,l/2;l,0;l/2,l/2 — 


Cl,l;l,l;0,0 1 
Cl,l;l/2,l/2;l/2,l/2 1 
Cl/2,l/2;l/2,-l/2;l,l "1 



TABLE III: Non zero structure constants of SU{2)\ and SU (2)2- 



by 



[Jq, 4>s,m\ = (S±m) 4>s,m±X- 



(9) 



With these commutation relations and those of the current modes (|6j), matrix elements of 
the form 

<<M<C ■ ■ ■ 4 a :0s', m (o,o)j\ ■ ■ ■ J\\M, (10) 

can be reduced to the structure constants, C Sl)mij]S2:m2]S3jm3 : 

2A S 



(0 S)S |</Vm(O)|</> s » s ») 



J ^ s,s;s' ,m;s" ,s" ■ 



(11) 



We list all non-zero structure constants, C Sl!mij . S2tm2 - S3>ma in Table III for SU "(2)i and SU (2) 2 . 

With the spectrum of the unperturbed SU(2)k model specified and the matrix elements 
of the perturbing field given, we are able to compute explicitly the matrix elements of the 
full Hamiltonian on the circle and represent it in matrix form. To be able to analyze this 
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Hamiltonian, we truncate (the truncation in the acronym TCSA) the Hilbert space by dis- 
carding all states with a chiral component whose level is greater than N tr . The resulting 
finite dimensional Hamiltonian matrix can then be diagonalized numerically, obtaining the 
spectrum of the perturbed theory. This procedure is particularly robust for relevant pertur- 
bations since the low energy eigenstates of the perturbed theory, |r), are localized on the low 
energy conformal states, {|c) Q } . Namely, expanding \r) into the conformal basis, {|c) a }, 



the coefficients b a , are primarily concentrated on the low energy conformal states \c) a de- 
termined by H su( 2) k . 

To extract physical quantities for the perturbed Hamiltonian Q, such as the mass gap, 
energy levels, correlation functions, etc., the choice of the system size, R, requires special 
consideration. For mR <C 1 (here m is the putative mass scale of the perturbed theory), 
the system lies in the UV limit where the conformal term, Hsu(2) k , of the full Hamiltonian 
dominates. In this regime the spectrum resembles that of the conformal Hsu(2) k where the 
energy levels scales as 1/R. In the IR regime, Rm ^> 1, the perturbation, J'dxQ(x), domi- 
nates and one expects a scaling of the form ~ i? 1 ~ 2As , where 2A S is the scaling dimension 
of the perturbing field. In general, the spectrum of the perturbed model must be extracted 
in a region of R where the conformal term and the perturbation are balanced in the sense 
that physical quantities remains stable under small variations of R. This region is usually 
denoted as the "physical window" . 

For theories where the dimension of the Hilbert space grows very fast, i.e. SU(2)k with 
k large, the truncation scheme proposed above may not yield accurate results. In those 
cases one can take recourse to a numerical renormalization group (NRG) improvement of 
the TCSA [2TH21)] . This procedure allows the TCSA to reach much higher truncation levels 
than that possible in its unadorned form. Taken together with an analytic renormalization 
group, it is possible to remove the effects of truncation altogether [21]. While this NRG 
has been tested extensively on relevant perturbations of conformal field theories, it has not 
hitherto been tried on marginal perturbations of CFTs. We will show in Section IV that 
the NRG can accurately predict the low lying spectrum even in the marginal case. 




(12) 



Q 
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III. SU(2) 1 PERTURBED BY THE SPIN-1/2 FIELD 



Our first test of the TCSA is the perturbation of the SU{2)\ WZW model by the singlet 
formed from the spin-1/2 operator of SU(2) 1} $ = (0i/ 2 ,i/2 4>i/2 ,-1/2 - 0i/2 -1/2 01/2,1/2): 

H rel = H su{2)l +h J dx . (13) 

The scaling dimension of $ is 2Ay 2 = 1/2, and so this is a relevant perturbation. Moreover 
the theory is invariant under h — > —h. This Hamiltonian is equivalent to the sine-Gordon 
model whose Lagrangian is given by [36] 



J d 2 a;(^(^) 2 + 2Acos(^)), (14) 



with f3 = 2tt. The correspondence between (13) and (14) is based on the identifications 



0i/2,i/2 (z) = exp P^^) ! 0i/2,-i/2 (z) = exp (^-i^^J , (15) 

which, up to phase redefinitions, are the bosonization formulas of the SU(2)i model [2]. 

We now demonstrate that the TCSA numerics reproduce the expected behavior of the 
sine-Gordon model at this value of the SG coupling. The spectrum of the sine-Gordon model 
at (3 2 = 2tt is composed of a soliton S and antisoliton 5* with mass M and two breathers Bi 
and B 2 with masses Mi = M and M 2 = a/3M respectively. The soliton, anti-soliton, and 
the first breather form a triplet under SU(2). The charges of the particles (S,Bi,S) are 
given by (1,0,-1) and they coincide with their S z quantum number. The second breather, 
B 2 , is a singlet under SU(2). 

Fig. [T] shows the low energy TCSA spectrum which reproduces the basic structure: a low 
lying triplet with mass M, and a single excitation at roughly a/3M. The expected value of 
the mass M can be determined from the coupling constant used in the TCSA. The relation 
between the coupling of the sine-Gordon model and the mass M is given as [37] 



h = A 



7lT(l 



SlT' 



2r (^) 



2-^ 



(16) 



where f = ^-^r = f . With h = 0.0942753.. we expect the mass to be 1 while from the 

1_ 8tt 

TCSA we find M = 1.0016, in excellent agreement. 
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FIG. 1: Plot of the TCSA data for the lowest six excited states in the S z = sector with the 
ground state energy subtracted. The lowest two excited states correspond to B\ (the S z = state 
of the triplet) and B 2 . The next four excited states are two particle states. The data are computed 
with a truncation level Nf r = 9. The fundamental triplet of particles can be seen to have mass 
M=1.0016. 

For the second breather, found roug hly at V3M, a more careful analysis of the TCSA 
data is required. For this excitation there are significant finite size corrections. These cor- 
rections can be understood as virtual processes on the cylinder [SI [18], which are suppressed 
exponentially as R becomes large. For B 2 these corrections are given by 

Am B2 (R) = -3V3e~ R/2 



d0 
2tt' 



e- mB i Rcoshe rn B2 Rcosh(9)(S*f t (9 + ztt/2) - 1) + 0{e~ CB ^ R ). (17) 



We can estimate the so-called error exponent to be a B2 = 1.105 [18]. Here S B f t {6) is the 
scattering S'-matrix between B 2 and one of the members of the triplet (S, B\, S) [38]. Fitting 
m B2 + Am B2 {R) to the TCSA data, we find that M B JM = 1.7322 ± 0.0003. Using this 
mass, we plot m B2 + Am B2 against the TCSA data in Fig. [2j which shows an excellent 
agreement between the theory and the TCSA data for R > 9. 

We now turn to the ground state energy E gs . The TCSA gives E gs as would be computed 
in conformal perturbation theory to all orders. This perturbative energy can be expressed 
as the sum of a linear term in R, proportional to a bulk energy density, ebuik, plus a term 

11 





mR 



FIG. 2: Single particle state B2: Analytic prediction for its mass with finite size effects (continuous 
line) compared to the TCSA data. 



given by the thermodynamic Bethe ansatz Etba [3 [39] . 

E gs = ebuik-R + E TBA (R). 



(18) 



The bulk contribution to E gs is given by [4TJ 



Cbulk 



M\ i 

tan - 

4 2 



-aM 2 R. 



(19) 



For sine-Gordon with f3 2 = 2tt, a = 0.14438.... The contribution from Etba{R) is given by 
the solution of a coupled set of integral equations involving the S-matrices of the various 
excitations in the model [5j [39] . At large R, this contribution reads 



/oo jn 
— cosh(#)e- 
■00 27T 



MRcosh(6») 



■V3M 



^ C0 sh(e)e-^ MRcosh ^ + O(e 

27T 



-2MR\ 



(20) 



and essentially marks the correction to the energy due to the spontaneous emission of a 
virtual particle from the vacuum which travels around the system before being reabsorbed. 
Fig. [3] shows the TCSA data against the theoretical values of E gs (including the full, not just 
the leading order large R, contribution coming from Etba{R))- There is a good agreement 
for R smaller than mR = 8 and then slight deviations thereafter, which can be reduced by 
increasing the value of (see Fig. |3j). 
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mR 



FIG. 3: Ground state energy: analytic prediction (continuous line) compared to the TCSA results. 
Black squares: TCSA with N tr = 9; orange circles: TCSA with N tr = 10. 

Above the single particle states, one encounters sets of two-particle states consisting of 
pairs of particles from the triplet and the singlet of the SG model. These two-particle 
states can be organized into SU(2) multiplets. For example, two-particle states involving 
the triplet decompose as (3 <g> 3) = (5 © 3 © 1). These states suffer finite size corrections 
due to scattering between the particles. These effects can be taken into account by solving 
the quantization conditions for the momentum in finite volume: 

2nrii = mRsinh.8i — i\n Sij(8i — 8j), i = 1,2, (21) 

where 0^ (i = 1,2) are the rapidities of the particles that parametrize their energy- 
momentum, (E,p) = (mcosh.6i,msmh.8i), and Sij is the scattering matrix between the 
two particles [H]. The solution of these equations, for a pair (m, ri2), yields the rapidities as 
a function of R, and so the energy of these states. Fig. [4] shows reasonably good agreement 
between the analytical and the TCSA results, particularly at large values of R. At smaller 
R, single particle virtual processes become important, leading to deviations between the 
TCSA and our analytic estimates. 
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2.8, 1 1 1 1 1 1 1 5 




FIG. 4: Two-particle states: comparison between the analytic results and the TCSA data. Left 
panel: a two particle state involving two triplet particles with total S z = and (m, 71%) = (— 1, 1). 
Right panel: a two-breather B 2 state with (711,77-2) = (—1,1). Continuous line: analytical results 



derived from Eqn. (21). Dots: TCSA data 



IV. SU {2)i PERTURBED BY CURRENT-CURRENT INTERACTIONS 

In this section we consider the perturbation of the SU{2)\ WZW model by the marginal 
current-current operator, 

H = H SU{2)1 -g Jdx J L {x) ■ J R {x). (22) 

Unlike the perturbation by the spin-1/2 field, here the sign of g matters as it differentiates 
the model's behavior in the IR limit [38] • For g < the perturbation is marginally relevant 
and asymptotically free. The corresponding theory is a massive integrable rational field 
theory (RFT) coinciding with the £77(2) Thirring model [12]. And for g > the coupling 
is marginally irrelevant and the theory undergoes a massless RG flow towards the SU (2)i 
fixed point. 

We address two questions here. We first ask if we can determine with the TCSA the 
universal correction (as explained below) to the ground state energy due to the marginal 
perturbation. And second, we question if the numerical renormalization group improvement 
of the TCSA works in the context of marginal perturbations [2TM25] . 
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FIG. 5: Plots of the ground state energy of SU(2) \ + Jl ■ Jr as a function of the marginal coupling 
g. Solid lines give the perturbative computation for both g > irrelevant (blue) and g < relevant 
(red). The black line shows the the second order perturbative correction in g. The data points 
represent the corresponding numerical data from the TCSA. 

A. Universal Term in Ground State Energy from the TCSA 



In the small g regime, the correction to the ground state energy on the cylinder, due to 
the marginal perturbation, is given in perturbation theory by 

.2 /o^\ 2 Q ^3 



-En 



6R C -¥ R [r 4 h -3l bR [R 



(23) 



where 



d 2 z\z - 1|~ 2 ; 



d 2 Zid 2 z 2 \zi — 11 ^ I ^ 2 — 1 1 l \ z \~ z 2\ 



(24) 



Our conventions are such that two point function of the perturbing current- current operator, 
<P= Jl - Jr, is 



(0(r)0(O)> 

while its three point function equals 

(0(ri)0(r 2 )0(r 3 )> = 



3 1 

4 \r\' 



kl2| 2 kl3| 2 |?~23| 2: 



(25) 
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Igl 

FIG. 6: Plots of the residual ground state energy (the TCSA data with the perturbative con- 
tributions, I 2 TCSA and I^ CSA , subtracted off) as a function of g for Nt T = 9 (red: relevant g; 
blue: irrelevant g). The average of the two results (black curve) is fitted with the function, 
g 3 (a + gbNtr{Ntr + 1)) (this average removes lower order logs that might appear from resumming 
high order contributions). The fitting parameters are a = 0.002 ± 0.001 and b = 94.12 ± 0.02. 

with b = 3/2. Although the theory is defined on a spacetime cylinder, under a conformal 
transformation the integrals can be written as integrals over the plane, as given above. 
The integrals I2 and J3 are UV divergent and require regulation. The divergent pieces of 
these integrals contribute to the non-universal piece of E (non-universal because their value 
depends on the regulation scheme). If e p is a short distance regulator in the plane (i.e. the 
difference of two plane integration variables is not permitted to be smaller than e p ), both I2 
and ^3 contain terms proportional to e p 2 . Under a conformal transformation, e p is related 
to a short distance regulator on the cylinder, e c , by e p = 2ne c / R. Thus the divergent pieces 
lead non-universal correction of the form 

R 

E^non-univ = («2fl' 2 + ^Q 1, ~\ ) (26) 

e 

c 

i.e. these corrections scale linearly with the system size. 

On the other hand, the universal corrections coming from I2 and J3 are independent 
of the details of the regulator. Those corrections are usually accessed through analytic 
continuation in the operator dimension x = 2, replacing the marginal perturbing operator 
with an operator with dimension x < 2. For sufficiently small x, the resulting perturbative 
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integrals become convergent. For such a range of x, one then introduces a regulator, e p , and 
expands in a Taylor series, finding that 7 2 and I 3 have the form [13J EH] : 

I 2 (x,e p ) = c 2 ef- a; )- 2 + / 2 (x,0) + O(e^); 

2— x 

I 3 (x, e p ) = c 3 ef-^- 2 + 6vr ^- J 2 (x, 0) + I 3 (x, 0) + O^). (27) 

Because everything is convergent, expressions can be developed for l 2 (x,0) and J 3 (x, 0) as 
a function of x. The universal terms at x — 2, are then analytically continued parts of the 
expansion of I 2 and J3 that are independent of the UV regulator (the non-universal terms, 
in contrast, are subsequently divergent when x is taken to 2). In the case at hand, the 
universal contributions are (131 El] 

I 2 (x = 2, e p ) univ = --(2 - x) = 0; 

e 2 "* 3tt 2 77r 2 

/ 3 (x = 2,e p ) univ = 6n-^I 2 (x,0) + I 3 (x,0) = -— - 2tt 2 = - — , (28) 

2 — a; 2 2 

where we see that 7 2 does not give a universal contribution, while J3 contributes from two 

terms in its Taylor expansion. 

An interesting question is whether this universal term can be extracted from the TCSA. 

While it has been suggested in Ref. (15] that the dependency of the TCSA data upon the UV 

regulator (here, the truncation level, N tr ) obscures such terms, it has been shown in Ref. [H] 

that upon the subtraction of leading and sub-leading order divergences, universal IR behavior 

can be observed. In order to extract such behavior, we need to compute the integrals 7 2 

and I3 with the TCSA regulator in place. This however can be straightforwardly done 

[15] . The essential idea is that these integrals can be expanded in powers of the integration 

variables Z\^- These expansions can then be compared with a Lehmann expansion of the 

corresponding correlation function. By truncating the Lehmann expansion to the same low 

energy states used in the TCSA, we are able to compute J 2i 3 with the TCSA regulator (see 

appendix A for details). The results are 

ll CSA = 3 -7rN tr (N tr + l); 

lJ CSA = 3n 2 N tr (N tr + l). (29) 

These expressions for J 2) 3 contain both non-universal terms (as is evident from the depen- 
dency on N tr ) and potential universal terms. The question becomes how to identify which is 
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which. We can do so by comparison with the evaluation of I2 and I3 in Refs. [HI HI], where it 
was argued straightforwardly that the contribution of I2 is purely non-universal. This corre- 
sponds to our finding that 7 2 is proportional to N tr (N tr + l) and indicates that the regulator, 
e p , used in Ref. |l3j HI] can be identified with the TCSA regulator via e~ 2 oc N tr (Nt r + 1). 
However lJ CSA is also proportional to N tr (N tr + 1), which seems to imply that the third 
order contribution to E gs , as with second order, is purely non-universal. 

To verify this non-universality of the ground state energy in the TCSA, we plot in Fig. 
[5] its value at small g for both signs, together with the perturbative results, obtaining ex- 
cellent agreement. To analyze more closely the possible presence of a universal term in 
the numerics, we plot in Fig. [6] the residual ground state energy arrived at by subtracting 
from the numerical data the perturbative contributions (solely non-universal) corresponding 
to I2 CSA + if" . To determine whether the numerics indicate any universal contribu- 
tion, we plot the residual as a function of g and fit the results to a function of the form 
g 3 (a + gbN tr (N tr + 1)). These fits put an approximate bound (the value of a) on the third or- 
der universal term consistent with the numerics. We find it to be considerably smaller than 
that found in Ref. [l3l HI] , consistent with our previous statement that in the regulation 
scheme used by the TCSA, there is no universal term at third order in the coupling. 

B. NRG and Marginal Perturbations 

In this subsection we apply the numerical renormalization group to the study of the 
marginal current-current perturbation of SU(2)i. The NRG is a technique that allows the 
TCSA to include states at much higher conformal levels than would be possible with a 
straight exact diagonalization. It does so by taking a cue from Kenneth Wilson's NRG [46J: 
it takes into account the states that have a weaker influence on the low energy eigenstates 
of the full theory, in this case high energy conformal states, only in numerically manageable 
chunks. It works in the case of a relevant perturbation because such perturbations guarantee 
that the high energy conformal Hilbert space only affects weakly the low energy sector of the 
theory. Thus it is not clear, a priori, whether the NRG will work in the case of a marginal 
perturbation where the high and low energy sectors of the theory are more tightly coupled. 
We will, however, see that the NRG does work, reproducing with high accuracy the results 
of the TCSA run with a straight exact diagonalization. 
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FIG. 7: Plots of the ground state (left), the first excited state (center), and the fifth excited state 
(right) energies as a function of RG step. The energies of the last RG step correspond to taking 
into account all states up to a truncation of Nt r = 11. Here g = 0.008 is marginally irrelevant. The 
NRG base matrix size and step size (N and A in the notation of [22J) are iV = 4767 and A = 500. 
The total number of states in the Verma module of the Identity is 14767. The solid line is the 
TCSA result done with an exact diagonalization at level Nt r = 11. 

Fig. [7] shows the RG evolution of the energies of the ground state, the first excited state, 
and the fifth excited state. The NRG results converge towards the exact diagonalization 
results with excellent accuracy. At least for the low lying energies in marginally perturbed 
conformal field theories, the NRG seems to be able to reproduce the expected energies, 
despite their dependencies upon the TCSA UV cutoff, N tI . 



V. SU(2) 2 PERTURBED BY THE SPIN-1 FIELD 

In this section we apply the TCSA to a deformation of the SU (2)2 WZW model [2j. The 
SU (2)2 WZW model is a c = | conformal field theory with three primary fields, the spin-0 
identity field, a spin | field, 0i/2,±i/2, and a spin 1 field, 0i,{±i,o}> with conformal weights 
A = 0, Ax/2 = jq, and A 1 = | respectively. 

SU (2)2 WZW can be thought of as a theory of three non-interacting massless Majorana 
fermions (the spin-1 fields). Perturbing SU{2) 2 with the spin-1 field, 

H = H SU (2) 2 + 9 j dx (0i,i 0i _i - 0i iO 0i,o + 0i,-i 0i,i) , (30) 

makes the fermions massive. Finding such a spectrum however is a strong check on the 
TCSA as the massive Majorana fermions are not simply expressible in the conformal current 
algebra basis. 
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FIG. 8: Plot of the energies of the four lowest excited states as a function of R in the sector S z = 0. 
The ground state energy has been subtracted. The data are found for a truncation level Nt r = 6. 
There is a single low lying state with mass M = 1.001 and spin S z = 0, which is a member of the 
fundamental triplet of particles. The higher energy level is three fold degenerate and corresponds 
to two-particle states. 

Fig. [8] shows the low lying excited states from the TCSA. The coupling g is fixed to 
a value where the mass of the first excited state is 1. The TCSA captures this state (the 
S z = state of the triplet of massive Majoranas) as well as a set of three two particle states 
with an energy approximately equal to 2. The Majorana fermions are non-interacting, hence 
their mass must scale linearly with g. This behavior is shown in Fig. [9^. We also consider 
the finite size corrections to one of the two-particle states. Even though non-interacting, 
the quantization condition for the two-particle state, Eqn. |2lJ is non-trivial as one needs to 
quantize with a doublet of distinct integers (ni, n 2 ), ni ^ n 2 . Fig. [9t> shows that the TCSA 



energy of the two-particle state matches the prediction derived from Eqn. 21 



VI. CONCLUSIONS AND DISCUSSIONS 



In this paper we have applied the TCSA to perturbations of SU(2)k- The general method- 
ology is illustrated with various examples. We began by studying the S77(2)i + Tr(g) per- 
turbation, equivalent to the sine-Gordon model at a particular value of its coupling. We 
showed the TCSA accurately captured both the ground state energy and excited state spec- 
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FIG. 9: a) Linear scaling of the mass of the triplet with the coupling g. b) A comparison of the 
energy of a two-particle state with (m,n2) = (—1,1) between TCSA numerics (black dots) and 
analytics (solid red line). 

trum including finite size corrections. We next turned to SU(2)i + Jl ■ Jr, a WZW model 
perturbed by a marginal interaction. We demonstrated that one can accurately identify the 
leading UV divergences that characterize the ground state energy computed by the TCSA. 
By subtracting these UV divergences (i.e. the non-universal contributions to the ground 
state energy) we isolated the universal contribution to the ground state energy due to the 
marginal perturbation. In the regularization scheme employed by the TCSA, we argued 
that the universal contribution at third order in the marginal coupling is zero. Finally we 
considered SU(2) 2 perturbed by Tr(g) 2 . This model is equivalent to three massive non- 
interacting Majorana fermions. The TCSA, is able to reproduce the expected spectrum. 
This is a non-trivial check of our methodology as the Majorana fermions do not have a 
simple representation in the current algebra basis employed by the approach. 

We have developed this capability to study perturbed WZW models in order to tackle a 
number of problems. In particular, we have two in mind. In the first, we plan to examine 
the dimerized-frustrated Ji — J 2 — 5 Heisenberg model whose Hamiltonian is 

/r = J>(i + <K-i) n ) s n ■ s n+1 + j 2 s n s n+2 . (31) 

Field theory analyses [27] argue that the low energy sector of this theory is equivalent to 

H = H sui2)l + g J dx J R - J l + h J dx (0i/ 2 ,i/2 ^1/2-1/2 - 01/2,-1/2 $1/2,1/2), (32) 

where g oc J 2 — J 2c and h oc 5. If the marginal perturbation g is absent, the model's spin 
gap, As, would scale simply with h: A s oc h 2 ^ 3 . In the presence of g, however, this scaling 
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is altered to become A 5 oc h 2 ^ 3 /\ logft)] 1 / 2 . However this altered scaling has been difficult 
to see in DMRG studies of the dimerized-frustrated Heisenberg model [28j EH]. It would be 
extremely interesting to analyze this scaling behavior using the TCSA. Our study here of 
577(2) 1 + Jl ■ Jr has thus lain the groundwork for this future study. 

The second problem that we intend to tackle is the study of possible integrable perturba- 
tions of SU(2)k for k > 1. There are indications, coming from Zamolodchikov-type counting 
arguments (17], that those perturbations do exists. We intend to study these perturbations 
with the TCSA, extracting both the spectrum of the model as well as evidence for or against 
their integrability. The example analyzed in this paper concerning SU (2)2 shows that this 
goal is within reach. 
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Appendix A: Evaluation of the Perturbative Integrals from the Current-Current 
Perturbation of SU(2) 1 

1. TCSA regularization and the Lehmann representation 

We implement the UV regularization explicit in the TCSA with a truncation of the 
Lehmann representation of the n-point functions |l5j 48J. As an explicit example, we com- 
pute the two point function using this truncation. Introducing projectors Vn tr an d PN tr that 
project out all chiral states of level higher than N tr , we can write the two point function 
appearing in the evaluation of the second order perturbative contribution to the ground 
state energy as 

(T#M)0(O,O)} = (TP Ntr P N ^x,t)P Ntr P Nt J(0,0)P Ntr P Ntr ), (Al) 
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where T is the usual time ordering operator. These projectors act to truncate the Lehmann 
mode expansion arising from the insertion of the resolution of the identity in between the 



fields of (Al): 



(0\<t>{x,t)\m) <8> \fh)(m\ g> (m|0(O, 0)|0). (A2) 

0<m,fh<N tr 

Here the truncated sum is over all states whose conformal level is less than N tr . Using the 
time and space translation operator, e ffi+,Pl , and defining z = e^* x+t ) and z = e^(~ tx+t \ 
we can rewrite the above as 

(O|0(O,O)(|m) g> \m)(m\ <g> (m\)(f)(0,0)\0}z m z™ } (A3) 

0<m,m<N tr 

so reducing the correlation function to a sum over powers in z and z. We see in this sum 
that no power greater than N tr of z or z appears. 

The case of the third order correction to the energy involving a three-point function 
follows the same strategy. We insert the projectors 

(T V Ntr V Ntr <t>(x, t)r Ntr V NtT (j>{0, 0)V Ntr V Ntr <f>{x', t')V Ntr P Ntr ), (A4) 

arriving at 

Y (O|0(x,t)|m) ® \fh)(m\ ® (fh\<f>(0,0)\n) ® |n)(n| ® (n\(j)(x',t')\0). (A5) 

0<m,»n,n,n<A r t r 

If we again use the space time translation operator and use the notation z = e^^ tx+t \ 
z = e^(~ tx+t \ w = e^( lx ' +t '\ and w = e^(~ %x ' +t '\ we obtain 

<0|^>(0, 0)|m> ® \fh) (m| <g> (m|0(0, 0) |n) <g> |n) (n| g> (n|0(0, 0) \U)z m z m 'w n w n , (A6) 

0<m,m,n,?i<A'i r 

that is a truncated polynomial at order Ntr i n -2, z,w,w. 



2. Evaluation of /2 



We want to compute the integral, ij , by regularizing as we did with (Al). Trans- 
forming the integral in I2 back to the cylinder we obtain, 



q^=2. { -- } j : : , (at) 



2.3 



where the factor of two counts the two contributions coming from time ordering. If we then 



expand the integrand as a power series in z, z and compare with Eqn. A3, we see we should 
truncate the sum as follows 

/27r\2 pR nO N tr Ntr 

/ 2 = 2.^-j J dx J dtJ2nz n ^mz™, (A8) 

whereupon the integral and sum are easily computed: 

/ 2 = 7riV tr (iV tr + 1). (A9) 

Comparing with the evaluation in [43, 44J, we see that we obtain a relationship between the 
TCSA regulator and the short distance cutoff, e p , used in [431 144] : 

e p = . 1 . (A10) 

y/N tr (N tr + 1) 

3. Evaluation of I3 

We now turn to computing lJ CSA : 

jTCSA = I dzdzdwdw -. -. -. — — -. (All) 

6 J {z - w){z - w){z - \){w - \){z - \){w - I) y ' 

Moving to the cylinder, the time ordering gives 6 contributions, t < t' < 0, t > t' > 0, 
t < < t' and those with t ■(-» t', all being equal. We show, as an example, the integral for 
t > > t'\ 

/2tt\ 4 f R , [ R , , [°° , f° ,. zzww , A „ n , 

UJ Vo " X /o " X 7o '"7 /'u- «-)U- „.)U- in,-- 1H,...M»-.. If (M2) 

where z = e^^ x+t \ z = e^^~ lx+t \ w = e^^ x ' +t '\ and w = e^^~ lx ' +t '\ We compute the 
contribution this makes by expanding the above in a power series of z, z, w, w and truncating 



it at the N tr -th power, matching the truncation in the Lehmann expansion in Eqn. |A6[ The 



integrand in (Al2) then becomes 





z 
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w 


2 






l 




l 




l 
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z — w 
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z 
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w — l 
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l — w/z 


2 


l - l/z 


2 


l — W 


2 


Z 


2 



(Al3) 



7 \ 7+1 / ? \ J+ 1 



j,k,lj,k,lEC 



V^ 1 z~' k - j -\ (A14) 

j,k,lj,k,l£C 
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where C is defined by 

< j + I + 1 < N tr ; 
0<j + k + l<N tr ; 
< J + r+1 < N tr ; 

< j + k + 1 < N tr . (A15) 
The integrals can now be done and the sum that so arises evaluated to be 

' ^(j+i+m + k + i)=* — 2 — • (A16) 

Once combined with the other five equal contributions, we obtain the final result: 

lJ CSA = 37T 2 N tr (N tr + l). (A17) 
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